Flux-Induced Vortex in Mesoscopic Superconducting Loops 
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We predict the existence of a quantum vortex for an unusual situation. We study the order 
parameter in doubly connected superconducting samples embedded in a uniform magnetic field. 
For samples with perfect cylindrical symmetry, the order parameter has been known for long and 
no vortices are present in the linear regime. However, if the sample is not symmetric, there exist 
ranges of the field for which the order parameter vanishes along a line, parallel to the field. In many 
respects, the behavior of this line is qualitatively different from that of the vortices encountered in 
type II superconductivity. For samples with mirror symmetry, this flux-induced vortex appears at 
the thin side for small fluxes and at the opposite side for large fluxes. We propose direct and indirect 
experimental methods which could test our predictions. 
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The question of where and under what conditions an 
individual vortex appears (or disappears) is still an active 
subjectLTQ. Here we consider an unusual setup: that of 
the Little-Parks experiment^ cl, i.e. a superconducting 
loop which encloses a magnetic flux. 

We find that, for a narrow range of enclosed fluxes in 
the vicinity of a half-integer number of quanta, a vortex 
is present in the loop. This vortex is completely differ- 
ent in nature from the traditional vortices known from 
type II superconductivity. In the traditional case, a min- 
imum local magnetic field is required for the formation of 
a vortex; here, its formation is governed by the enclosed 
flux. As the magnetic field is increased, instead of hav- 
ing more and more vortices coming in, the present vortex 
disappears (and a new vortex appears only after an addi- 
tional quantum of flux is enclosed by the loop) . Whereas 
traditional vortices cannot exist along films which are 
thinner than.-the coherence length (they "need room for 
their cores" )□, the present vortex exists for arbitrarily 
thin shells. Another atypical feature of the present vor- 
tex is its high anisotropy, which depends on the thickness 
of the shell. Moreover, there are cases in which the order 
parameter vanishes on a surface rather than a line. 

Additional differences arise from the assumption that 
the total magnetic field remains equal to the applied field: 
whereas traditional vortices exist for type II supercon- 
ductors, for which the Ginzburg-Landau parameter n is 
larger than 2~ 1 / 2 , the present vortex is insensitive to the 
London penetration depth and bears no connection with 
k; whereas traditional vortices cannot stay in eouilib- 
rium if the distance to the boundary is too shortu, the 
equilibrium position of the present vortex is a continu- 
ous function of the flux and does reach the boundaries. 
This assumption is justified when either the thickness or 
the width of the loop is small compared to the magnetic 
penetration depth. The case in which the induced field 
is not negligible will be considered elsewhere. 

We shall consider two complementary geometries: 
samples which almost have axial symmetry, which will 
be analyzed using perturbation, and samples with ec- 
centric cylindrical boundaries, which will be studied by 
means of variation. In the first method the zeroth or- 
der configuration is a cylindrical shell of superconducting 



material, with inner radius R, outer radius R' = R + w 
and heigth h, where R, w and h are constants. In early 
experiments!!] h was,much larger than R and w, whereas 
recent experimentstlt] meet the opposite situation. We 
consider a uniform magnetic field parallel to the axis of 
the shell and study the onset of superconductivity. In 
this regime the magnetic field is still unaffected by the 
supercurrents and the Ginzburg-Landau equation is lin- 
earized to the eigenvalue problem 

Ho^Po = (tV - 2vrA/$ ) 2 O = W£o = Eo^o- (1) 

■00 is the order parameter, A is the magnetic poten- 
tial, $o is the quantum of fluxoid (<!>o > 0) and £o is 
the coherence length at the onset of superconductivity. 
The notations for the function ipQ, the operator Hq and 
the eigenvalue E have been introduced to remind of the 
eigenf unction, the Hamiltonian and the energy in a per- 
turbation problem. Eq. ([!]) is subject to the condition 
that («V — 27rA/$o)V' has no normal component at the 
boundaries. (For the moment, tp = ipo.) For a cylindri- 
cal shell and cylindrical coordinates (r,6,z), Eq. (|l|) is 
separable and has been solved in [ fill . Writing i/jq = 
TZ(r)Q(0)Z(z), one easily obtains Z[z) = cos(kTrz/h) 
and <d(9) — e~ mie , with k and m integers. Since Eq 
always increases with fc 2 , we pick k = and the solution 
of Eq. (||) reduces to 
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The winding number m is chosen so that lowest value of 
Eq is obtained. For a thin shell it is the closest integer to 
the number of flux quanta enclosed by the shell. lZ m was 
obtained in [ it depends on the magnetic field and 
involves a combination of Kummer functions, adjusted to 
meet the boundary conditions. 

A salient feature of the solution (S) is that there exist 
magnetic fields for which Eq. ([!]) is degenerate. In the 
limit of a thin shell this occurs when the enclosed flux $ 



equals 



^)<E>o : for this flux, the eigenvalues Eq ob- 



tained for m and for m + 1 are the same and the lowest 
possible. If the shell is not thin, we can still define $ as 
the magnetic flux enclosed by the "representative circle" 
r = R = R + w/2. There still exist fluxes (usu- 
ally close to (m + |)$o) where the values of E for O = 
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7£m (?") exp(— miff) and for ipo = 1Z m +i(r) exp(— [m+l]i9) 
coalesce. As the magnetic field is swept accross the situa- 
tion $ = (f>*„, the order parameter ipo changes discontin- 
uously. With it, there is a jump in measurable quantities 
such as the current around the shell. 

We perturb now the problem and consider a sample 
which is not a perfect cylindrical shell. To fix ideas, we 
still keep h and R constant, but take a nonuniform width 



D(9) 




(3) 



where w is the average width, e C 1 and /?_,■ = j3j . Other 
deviations from the perfect cylindrical problem lead to 
similar behavior. 

The "Hamiltonian" is still Ho, but the eigenvalue and 
the eigenfunction are perturbed by the change in the 
boundary geometry. We write ip — ipo + eipi + e 2 ip2 + ■ • • 
and E = Eq + tE\ + e 2 E 2 + ■ ■ ■, where ip is the order pa- 
rameter and E~ x / 2 the coherence length. The eigenvalue 
in our problem is related to the temperature T by 



R 2 E = (T C -T)/(T C -T R ), 



(4) 



where T c is the critical temperature in the absence of 
magnetic field and Tr the temperature at which the co- 
herence length equals the internal radius R. The only 
reason for the requirement of a "mesoscopic" sample is 
to decrease Tr and thus widen the temperature scale 
T c — Tji] our formalism is valid for arbitrary positive R, 
w and h, and may cover the entire range from a very thin 
ring to an almost full disk. In order to decrease Tr we 
also require clean materials (usually type I), which have 
a large coherence length at zero temperature. 

To proceed, we define a metric by integrating in the 

unperturbed region: (fa, fa) = fodBfZ rdrfafa. It fol- 
lows that 



i,H <p 2 ) - (Ho4>\,4> 2 ) = R' J f^ 2 ~^~T ~~ ^"^T 



d6. 



(5) 



Far from we obtain the sequence of equations 
(H - E )tpi = Eiip , (Ho - E )ip 2 = Eii>i + E 2 ip , etc. 
The difference between this procedure and standard per- 
turbation theory is that Ho is not Hermitian; instead, it 
obeys Eq. (|J). To take this into account, we need dipi/dr, 
evaluated at r = R' . To obtain this we expand ip around 
r = R' and require to all orders of e that the normal 
current vanishes at r = R + D(0). In this way we obtain 
that Ei — 0; the expressions for tpx and E 2 are lengthy 
and will be reported elsewhere. 

For $ w the direct perturbation scheme described 
in the previous paragraph would lead to a divergent ip\. 
This divergence is due to the degeneracy at $ = <&* m . In 
this case we use degenerate perturbation theory, i.e. we 
write 



ipo = TZm(r) exp(-miO) - jH m+ i(r) exp(-[m + l]i0), 



with the normalization lZ m (R) = TZ m +i(R) = 1 and 7 
a coefficient which still has to be determined. Writing 
2irR 2 A = $ [6 m + eS]r§, with R 2 b m = $^/$o, and sub- 
stituting into (|J) <pi by ipo and <p2 by TZ m (r) exp(— miff) 
or lZ m+ i(r) exp(— [m + l]i6), leads to a system of equa- 
tions for 7 and E\. 



A m S + B m E\ = C m ^i7 , 
A' m+1 S + B m + X E X = C m hll , 



where 



R' 



A m =2 (mr- b m r s )K 2 m dr , 

JR 
f R' 

A' m = 2 (mr- b m ^r 3 )TZ 2 m dr , 

JR 



B„ 



rlZ 2 dr 



(7) 

(8) 



(9) 
(10) 

(11) 



(6) 



C m = \wR"R m (R')n m+l (R') x 

[T - (b m R' - (m + l)/R')(b m R' - m/R')] (12) 

and R has been taken as the unit of length. This system 
has two solutions and we choose the one with lower E\. 
Note that this time E\ ^ 0. After obtaining ipo and E±, 
ipi and T 2 were also evaluated, but the details will be 
reported elsewhere. 

The value of E provides the temperature at which the 
shell becomes superconducting. We have calculated it for 
a sample with a reasonable experimental shape: D(9) = 
0.2R - O.OLRcos(0). The results are shown in Fig. 0. 
A customary approximation (e.g. [ ^)) for results in a 
thick shell is obtained by averaging the one-dimensional 
result for ideally thin shells. This procedure gives T av = 
J (nBr/^o — m/r) 2 dr/w, where B is the magnetic field, 
m is the integer that minimizes this expression and the 
integral is from R to R' . A somewhat better approxima- 
tion is obtained by taking lZ m constant (e.g. [||]). This 
gives E const = f(irBr/$ Q - m/r) 2 rdr/ J rdr. T const is 
also shown in Fig. [j] for comparison. We see that T const 
is quite a good approximation, but it exhibits cusps when 
the value of m changes, whereas for our treatment T is 
smooth. This qualitative difference does not arise from 
the dependence of lZ m on r; what happens is that the 
superposition of two values of m produces a sinusoidal 
dependence of \ip\ 2 on The size and the extent of the 
decline of T c — T near $^ are proportional to e(3\ , which 
measures the deviation from a uniform shell. 

In Fig. ^ we show contour plots of \ip\ in a narrow radial 
strip next to 9 — for six equidistant fluxes between 
$ = 0.501364$ and $ = 0.501382$ - % = 0.501376<J> 
is contained in this range. We have chosen an example 
with the width D symmetric about 9 = and 9 = ir, so 
that the plots can be continued by placing a mirror at 9 = 
0. As $ increases and comes close to $qi IV'I decreases 
at the axial line (r = R + D(0), 9 = 0), until in frame (b) 
ip(R + D(0), 0) = 0: this is the field at which the vortex 
appears. As the flux is increased further, the equilibrium 
position of the vortex moves towards the inner boundary, 
which is reached in frame (e). If the field is increased 
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further, the vortex disappears. The extent of the flux 
range for which the vortex exists is proportional to e(3±. 

In the general case the vortex is located at 9 = 
arg(— C m /?i); if m is small and D is symmetric, as in the 
example we picked, this corresponds to the thinnest part 
of the shell. For large values of m (m > 5 for w = 0.2R), 
C m becomes negative and the vortex is located at the 
opposite side of the shell. 

Knowing ip(r,6), we can evaluate the supercurrent. 
Fig. H shows the streamlines for the same fluxes consid- 
ered in Fig. ||. Each frame shows 1/20 of the shell, with 
the plane 6 = at the left. 

The perturbational approach is theoretically instruc- 
tive, but since the range of existence of the vortices in- 
creases with nonuniformity, we are interested in sam- 
ple shapes that are very far from axial symmetry. We 
conjecture that the scenario found above is generic for 
samples that are not symmetric. To support this con- 
jecture, we consider a sample with cylindrical bound- 
aries. The inner (outer) radius will again be denoted 
by R (R 1 ), but the axes of both boundaries do not coa- 
lesce, will still be the flux through a circle of radius 
R — ^(R + R 1 ). It is natural to use bipolar coordinates, t3 
x + iy = ctanh[(a + ij3)/2], where x and y are the 
Cartesian coordinates perpendicular to the axes and c 
a constant which defines the lengthscale. The lines of 
constant a are circles; we assign a\ = arcsinh(c/i?'), 
a-2 = arcsinh(c/i?), and the sample occupies the re- 
gion ai <a<a2, — it < (3 < n. Thinner samples 
(smaller (R' — R)/R) are obtained by taking small dif- 
ferences a2 — ai, whereas more eccentric samples (larger 
D(w)/D(0)) are obtained by taking both a\ and «2 small. 
The shape shown in Fig. ^ corresponds to a.\ = 0.5, 
a 2 = 0.7. 

The magnetic potential can be written in the form 



$(cosh a + cos 0) 



f(a,P)-f(a,P) 



R 2 



/(<*,/?) 



ircR 2 

c 2 sin -2 /? [sinh a/(cosh a + cos (3) 
2 arctan(tan ^ tanh ^) cot /?] , 



0, 



(13) 



with a — arcsinh(c/-R). 

An approximation for the order parameter ip -sszjill be 
obtained by means of a variational procedure.Ej The 
eigenvalue and the eigenvector for the onset of super- 
conductivity may be obtained by minimizing the ratio 
JT |(*V — 2nA/<f> )4>\ 2 dV/ f \ip\ 2 dV, where dV is the el- 
ement of volume and the integrals are over the sample. 
Using the boundary conditions, ip can be written as a 
Fourier series 



oo oo 



ip = ^2 ^2 ( ° mi cos ( liT 



m= — oo 1=0 



2a — ct2 — a\ 
a 2 — ai 



s m i sm 



(21 + 1)tt 2a - a 2 - ai 



a 2 — ai 



(14) 



("0 is again independent of z.) We truncate now this 
series and leave only a finite set of nonzero coefficients. 
The number of terms which is required to achieve some 



desired accuracy increases with $ and with (R' — R)/R, 
while it is not very sensitive to eccentricity. In our calcu- 
lations we have kept eight coefficients: c mo _2,o, c mo _i,o, 

Cm o ,0, C mo + i : o, Cm +2,0, Sm -1,0, Sm o ,0, Sm o +l,0j where 

the integer mo is chosen to obtain the minimum eigen- 
value. 

Regarding these coefficients as a vector 

Cmo+2,0, Smo-1,0) S mo ,0: Wl,o}i 

we can write /|(iV - 27rA/$ )?/'| 2 e?F = v'Mv and 
J \tfj\ 2 dV — v^Qv, where the elements of the matrices M 
and Q are easily evaluated. We are therefore left with 
the minimization of the ratio (v'Mv)/(v'Qv). This ra- 
tio is minimized for the eigenvector of Q -1 M that has 
the lowest eigenvalue. The error introduced by trunca- 
tion may be estimated by the size of the discontinuity of 
-0 at the values of $ for which mo changes. 

Using this method, we rediscover the behavior ob- 
tained for almost axially symmetric samples: near half- 
integer fluxes, a vortex appears at the outer boundary 
and moves towards the inner boundary as the flux is in- 
creased. What looks more striking, is that again these 
vortices enter through the thin side for small numbers of 
quanta and through the opposite side when <&/ <&o is suf- 
ficiently large. If are chosen to reproduce the shape 
used in our previous perturbation example, the length of 
the range for which the vortex is present agrees to about 
10~ 8 <I>o with the result obtained by perturbation. The 
values of the flux for which the vortex appears agree to 
about 10 _6 <E > o m both methods. Fig. [| shows contour 
plots of |^| for a-y = 0.5, a2 — 0.7 and for fluxes chosen 
so that the vortex is approximately midway between the 
inner and the outer boundary. For these values of u\ t 2 
and for $ < 5$o> we estimate the error in tp to be a few 
percent of the maximum value of \ip\. 

In conclusion, we have found a setup for which a new 
kind of vortex exists. Its existence is governed by the 
flux enclosed in the entire shell rather than by the local 
magnetic field. Likewise, the position of the vortex is a 
function of the entire sample (and of the flux) rather than 
of local defects. The ocurrence of this vortex is a nearly 
periodic function of the flux. The presence of a vortex 
implies that the flux lies within a very narrow known 
range. 

It should be stated, however, that the magnetic field in 
the superconducting material itself does influence the na- 
ture of these vortices. For instance, we have found that if 
the field is constant in the region < r < R and has op- 
posite sign outside this region, then, as $ increases close 
to $* n , a vortex and an antivortex form at a point some- 
where in the middle of the shell; as the flux increases fur- 
ther, the vortex moves from this point towards the inner 
boundary and the antivortex towards the outer bound- 
ary, until they disappear. Moreover, if the field vanishes 
in the entire shell, then the vortex becomes a cutf3 i.e. 
the order parameter vanishes on a surface that connects 
the inner and the outer boundaries. 
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Finally, we discuss the experimental possibilities for 
the detection of the vortex predicted here. We require a 
sample such that its size is of the order of the coherence 
length for a reasonable range of temperatures. Typical 
experiments for this purpose use samples made of Al, 
with a perimeter of a few /im. In order to observe the 
vortices directly, we require an imaging technique with 
spatial resolution of the order of 0.1 /im. This xequire- 
ment is met by japanning tunneling microsccipy,ll3 mag- 
netic decorationEj and electron holography.EZl When in- 
terpreting electron holography it should be born in mind 
that for the vortex predicted here only the fluxoid equals 
its magnetic flux alone will be much smaller and tem- 
perature dependent. Clearly, in order to sense the local 
density of states or the current density, the size of the or- 
der parameter must be above some threshold required by 
the sensitivity of the experiment. In this regime Eq. (|l|) 
for the onset of superconductivity becomes insufficient 
and we have to analyze the nonlinear Ginzburg-Landau 
equations. We have partially analyzed this situation and 
anticipate here some results. The existence and behavior 
of the vortex remain qualitatively the same as reported 
here for some temperature range below the onset of su- 
perconductivity. This temperature range increases with 
the deviation from a uniform shell. 

Direct observation of the vortex is expected to be a 
difficult experiment, since the requirement of large co- 
herence lengths implies a small gap. In addition, evap- 
orated Al usually has a rough surface. Therefore, it 
might be helpful to obtain indirect evidence for the sce- 
nario encountered here. In the limit of thin shells, we 
know that there exist critical points at fluxes close to 
and at some temperature which we denote T 2 M . Be- 
low T 2 (m) the current / around the loop exhibits hystere- 
sis, whereas above T^ 1 it is reversible and continuous. 
Moreover, above and near T 2 , dl/dQ has a maximum 
near $ = which diverges as is approached. EH 

Our analysis of the nonlinear Ginzburg-Landau equations 
shows that these features remain valid when the shell is 
thick. In particular, the divergence of dl/d$> at the crit- 
ical point is not smeared. 

We therefore need a technique that measures dl/d&, 
which is proportional to the ac magnetic susceptibility. 
This is measured by applying a bias flux with a super- 
imposed small ac signal. The response is the ac compo- 
nent of the induced magnetic flux in some region. This 
induced flux has already beep, measured by means of 
a SQUID microsusceptometerB Other techniques which 
seem appropriate for this measurement are the piezoresife 
tive cantilevertJ and the ballistic Hall magnetometer .E2l 
The Hall magnetometer seems particularly promising, 
because, if only part of the sample is located on the ac- 
tive region of the probe, it would allow to measure the 
flux induced at that part. 
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FIG. 1. Temperature T for the onset of superconduc- 
tivity as a function of the magnetic flux <E> enclosed by the 
"representative" circle. The graph shows two curves, which 
nearly coincide. The lower curve corresponds to the formal- 
ism developed here; for comparison, -Econst (see text) has also 
been drawn. The inset shows an enlarged view of the region 
enclosed by the rectangle. The lower curve in the inset looks 
smooth, but was calculated using three different algorithms in 
different regions: m = 0, degenerate perturbation and m = 1; 
likewise, near $ = $1 w 1.5<£>o, m = 1, degenerate perturba- 
tion and m = 2 match smoothly, and so on. 



d e f 

FIG. 2. Contour plot of the absolute value of the order pa- 
rameter in a narrow radial slice of the superconducting shell, 
for increasing values of the magnetic field. The horizontal co- 
ordinate is \0\, in the region < \9\ < 0.00027T; the vertical 
coordinate is R < r < 1.19i2. The plane 6 = is located 
at the left of each frame. In (a), has a minimum at the 
line (r = 1.19-R, 6 = 0). In (b), ip vanishes along this line in 
the outer boundary and the vortex appears. As the magnetic 
field increases, the vortex shifts towards the inner boundary 
r = R and reaches it in (e). In (f), \tp\ > everywhere and 
the vortex has disappeared. 
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FIG. 3. Streamlines in a piece of the superconducting 
shell, for the same magnetic fields as in Fig. |^. The coor- 
dinates are the same as in Fig. but now < \6\ < O.lir. 
In (a) and (b) we see streamlines through 9 = that carry 
electric current to the right (assuming that the magnetic field 
comes out of the page). These streamlines circulate around 
the superconducting loop (enclose the hollow region). At the 
right we see part of the screening currents, in clockwise cir- 
cuits that do not enclose the hollow region. In (c) we see 
two counterclockwise streamlines around the quantum vortex, 
three streamlines around the loop and the screening currents 
at the right. Along a path around the vortex the phase of 
the order parameter changes by 2n. In (d) the vortex has 
moved towards the inner boundary; the current around the 
loop has now become counterclockwise. As the magnetic field 
increases, the vortex moves down until it disappears and the 
counterclockwise current around the loop increases. The bot- 
tom frame shows a "panoramic" view of half of the shell, 
< \0\ < it. The magnetic field is the same as in (d). The 
thick lines are streamlines and the thin lines are lines of con- 
stant 1-01. On this scale, |V>| looks independent of r and the 
currents around the vortex and around the loop do not show 
up. 



FIG. 4. Contour plot of \ip\ for a sample with cylindrical 
inner and outer boundary. Darker areas denote smaller values 
of \ip\. For the chosen flux values, ip = roughly in the middle 
of the darkest areas, (a) $ = 2.5713$ - (b) $ = 3.61$ . The 
flux range for the presence of the vortex in case (b) is about 
forty times larger than in case (a). 
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